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I. INTRODUCTION 

Magnetic fields are present in most astrophysical sys- 
tems [Hi! |a 11 1* but their origin is still unknown. Spiral 
galaxies typically contain magnetic fields of about 1CP 6 G 
that are aligned with the spiral density waves 0. A plau- 
sible explanation is that these fields have been produced 
from the exponential amplification of an initially weak 
field by a mean- field dynamo [fj, [?J , in which a seed-field 
was amplified by the differential rotation of the galaxy 
in conjunction with magnetohydrodynamic turbulence. 
Magnetic fields of comparable amplitude are also found 
on cluster scales ||. 

The seed-field strength required at the time of com- 
pleted galaxy formation for a galactic dynamo to produce 
the present magnetic field strength of a few aG is usually 
quoted in the range - (1CT 23 - 1CT 19 ) G [f], H H S Hf 
These estimates, though, are obtained by considering the 
dynamo amplification in a flat Universe with zero cos- 
mological constant for "typical" values of the parameters 
of the aw-dynamo. According to Ref. [[| these lower 
bounds can be relaxed to about 1CU 30 G, for a Uni- 
verse with a dark energy component (e.g. a cosmological 
constant or quintessence), which appears to be favored 
by recent results from high-redshift Supernovae obser- 
vations, cosmic microwave background experiments and 
large-scale structure data. 

Most proposed models of primordial seed-field genera- 
tion either fail to meet these requirements or invoke non- 
standard coupling between the electromagnetic field and 
the inflaton scalar sector responsible for a primordial pe- 
riod of accelerated Universe expansion (see, e.g. Ref. Q 
and references therein). 

The most economic and conservative physical mech- 
anism for producing such a seed- field was proposed by 
Harrison [l(J]. The mechanism relies on the fact that 
weak magnetic fields are generated during the radiation 
era in regions that have non-vanishing vorticity. The 
main problem with Harrison's mechanism is that the re- 
quired vorticity of the plasma does not have any dynam- 



ical origin, but has to be put in by hand, as an initial 
condition. In this paper we show that rotational veloc- 
ity of the plasma is unavoidably induced gravitationally 
by the nonlinear mode-mode coupling of primordial den- 
sity perturbations generated during inflation, arising at 
second order in perturbation theory. Therefore, we pre- 
dict that a minimal seed magnetic field is generated by 
the differential rotational velocity of ions and electrons. 
The resulting magnetic field power-spectrum is fully de- 
termined by the amplitude and the spectral index ns of 
primordial scalar perturbations. The rms amplitude of 
the seed-field at recombination is B 10 _23 (A/Mpc) -2 
G, on comoving scales A > Mpc. 

The plan of the paper is as follows. In Section II we 
briefly recall Harrison's mechanism and we show how the 
seed magnetic field is produced. In Section III we com- 
pute the power-spectrum of such a seed-field in terms 
of the power-spectrum of primordial scalar perturbations 
generated during inflation. Finally, in Section IV we draw 
our conclusions. 



II. GENERATION OF SEED MAGNETIC FIELD 

Let us first recall how Harrison's mechanism 0] works. 
Consider a rotating region in the expanding early Uni- 
verse (with scale factor a), consisting of matter (mostly 
protons) with average energy density p m , electrons and 
photons with energy density p 7 . Electrons and photons 
are tightly coupled and are considered as a single fluid 
m Ref. 101. Let uj m and tt> 7 be their angular velocities. 
In absence of interactions, their angular momenta, pro- 
portional respectively to p m uj m a 5 , and p 7 w 7 a 5 , are sepa- 
rately conserved (see also Appendix 0. As the Universe 
expands, they scale differently, uj m tx a~ 2 and oj 1 oc a . 
Ions and the electron-photon fluid spin down at different 
rates and then currents and magnetic fields arc gener- 
ated. 

A formal derivation of the resulting magnetic field 
strength in terms of the vorticity in the proton fluid in- 
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eluding the electron-proton coupling can be found in . 
It is necessary to consider the evolution of the multi- 
component fluid formed by photons, electrons and pro- 
tons during the radiation era. At temperatures T > m e 
(where m e is the electron mass) , the interactions between 
ions and electrons are strong due to copious lepton-pair 
production, and they are locked together. The interac- 
tion between electrons and photons is also very strong. 
This means that all the plasma has the same angular 
velocity at T m e and no magnetic field can be gener- 
ated. Below this temperature, electrons and photons are 
still tightly coupled through Thomson scattering. Pro- 
tons and electrons are tightly coupled through Coulomb 
scattering (scattering between photons and protons can 
be neglected as the coupling is weaker) and so the pho- 
ton fluid drags the protons in its motion. The differ- 
ence in the strength of these interactions and in the mass 
of electrons and protons leads however to a small differ- 
ence in the electron and proton fluid rotational velocities, 
that gives rise to non-vanishing currents and magnetic 
fields. In our case the rotational velocity of the plasma 
is sourced by the non-linear evolution of scalar (density) 
perturbations. The generation of the magnetic field ends 
at recombination, when electrons and protons combine to 
form neutral hydrogen and radiation decouples from mat- 
ter. This means that the generation of the magnetic field 
starts around T ps m e and ends when T ss T rec . Dark 
energy and neutrinos can be ignored throughout the evo- 
lution. Moreover, during the radiation era we can safely 
neglect the role of the Cold Dark Matter (CDM) compo- 
nent, that is, we consider a plasma with only baryons 
(protons), electrons and photons, and electromagnetic 
(EM) fields. After matter-radiation equality, CDM plays 
the dominant gravitational role and it enters our equa- 
tions determining the evolution of perturbations. 

Let us now go into the details of this mechanism. The 
momentum equation for the interacting components can 
be written using the total energy-momentum conserva- 
tion equation — The energy- momentum tensor of 
each component is not conserved independently and its 
divergence has a source term that takes into account the 
energy and momentum transfer among the components, 
T$ v = Qi A \ with A = f,p,e,EM. We will describe 
the proton, electron and photon components as approxi- 
mately perfect fluids, thus their energy momentum tensor 
can be written as 

= (p( Q ) + p^)u^u {a)v + p^Si! , (1) 

with a = e,p, 7 and p^ = w^ a > p( a \ We can expand the 
energy density and four velocity of each component as: 

u (q) " = - (S% . (2) 

In what follows we will neglect the electron and proton 
pressure, i.e. we set ~ « 0, while for the pho- 
tons we have u/ 7 ' = 1/3. 



The crucial step in our computation is that we have to 
include second-order terms in the metric and the matter 
perturbations 0, 0, ^| , since rotational component of 
the velocity, and hence magnetic field, are only generated 
at second order. Indeed linear vector modes are not gen- 
erated during inflation and, by Kelvin's circulation theo- 
rem vorticity is conserved for a perfect fluid. Of course, 
this conservation holds at second order as well, and it 
applies to the vorticity of the fluid which generates the 
gravitational field. As explained in detail in AppcndixlAl 
this conservation law does not prevent the occurrence of 
a rotational velocity component of the charged particles, 
thus giving rise to a non-zero magnetic field. 

Let us expand = 6$ + 5$ and " = vjg 11 + 

u (2) ' wnere ^ ne superscripts (1) and (2) indicate first 
and second-order perturbation quantities. We will work 
in the Poisson gauge, which is the generalization of the 
so-called longitudinal (or Newtonian) gauge when vector 
and tensor modes are allowed. The perturbed metric 
reads 

ds 2 = a 2 (ri) {-(1 + 2cf>)dri 2 + 2 Xt dr]dx l + 

+ [{l-2ip)8 ij + X ij]dx i dxi}. (3) 

In this gauge, the first-order scalar perturbations are 
equal if there are no anisotropic stresses, <jP-> — ijM' = p. 
The vector perturbation is Xii an d we will assume that 
it is not generated at first order x^p — 0- We will also 
neglect primordial tensor modes. 

Although first order primordial vector and tensor 
modes are absent, the non-linear evolution of the pri- 
mordial scalar perturbations gives rise to non-vanishing 

vector and tensor perturbations (and corrections to the 

(2) 

scalar perturbations) at second order, Xi = Xl an d 

Let us consider the momentum continuity equation for 
each fluid component (whose general expression up to 
second order in perturbation theory is reported in Ap- 
pendix [SJ. The contribution of the electromagnetic field 
can be included in the source of the charged fluid compo- 
nents, as Tl% M)v = FrfjP, with f = en(uW - u^ 13 ), 
where we have assumed charge neutrality (n — ~ 
72^). The momentum equation for photons up to sec- 
ond order can be written as 

^f- (K (7) + X*)' + \di5^ + M + - 

_ + 1(^(7)2) . + 1 ^ = K? | (4) 

and for protons 

P { o P) ((v^+Xiy + Hiv^+xA + d^ 

- 2 ^ + i(^) 2 ) >i + i(Ai) 

= en(E i + e tJk vWB k )-K?, (5) 
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where TL — a' /a is the expansion rate in conformal time. 
The momentum equation for electrons is similar to that 
for protons, but it has the opposite sign in the righ-hand 
side terms and an additional source, — k^ 7 , due to the mo- 
mentum transfer between the electron and photon fluids. 
This momentum transfer is due to Thomson scattering 
and is given by n e f — — | pg 7 V'(«j 7 — v- ), where the dif- 
ferential optical depth is r' = anax, and gt is the Thom- 
son cross section. The momentum transfer between elec- 
trons and protons is due to Coulomb scattering and can 
be written as k^ p — —n(v^ — v\ e ^)/T e , where the colli- 
sion time between electrons and protons is r e = m e a/ne 2 
in terms of the conductivity of the plasma a. 

During the whole period we are interested in, up to 
the time of hydrogen recombination, the collision time 
between electrons and photons (oc r' _1 ) and the electron- 
proton one are much shorter than the Universe expansion 
time, thus momentum transfer is very efficient. A tight 
coupling expansion of the momentum equations gives, to 



,(7) 



,0) 



We can obtain an equation for the vorticity of the pro- 
ton fluid by taking the curl of Eq. JSJ) and combining it 
with Maxwell's equations 



d_ 

drj 



(p) 



fli H B 



e 
in. 



= 2e ijk a^v^ k +- 



(6) 

where we have kept terms up to second order in the met- 



ric perturbation and defined £ 5 



(p) 



jk d j V^ k /a and 



Qi = Cijkd 3 x I O" I n the last equation the diffusion term 
can be dropped in the highly conductive protogalactic 
medium. 

The seed magnetic field can be written as 



B, 



e 



,(p) 



dfjae ljk d^'v^ k 



(7) 

where the subscript "7" denotes the initial time, corre- 
sponding to T ~ m e , when Bj = 0, and where we intro- 

(A) (A) 

duced the useful combination p\ = Q + fij. 

Taking now the curl of the photon momentum equation 
ijljl we obtain an equation for the photon vorticity 



M 7) )' = Jew* 



9 3 y 



(8) 



where we have set = v. Solving Eq. JHJl we finally 
get: 



= ^«' + ^ J (-8* <p&&™ 

_ ~v j d k d l v l + 4 ^ip'v k \ (9) 
where we have set (3^ ~ (3^ = (3. 



Equation J7J with the solution Eq. © is our main re- 
sult and provides the resulting magnetic field in terms of 
the metric and matter perturbations. In order to finalize 
our computation we have to consider the linearly per- 
turbed Einstein's equations. For a Universe dominated 
by a fluid with equation of state p = wp, with w = const., 
we can express the fluid velocity t> W 1 and the energy 
density perturbations #W in terms of the gravitational 
potential tp as 



<5« = 



3(1 + w)H 2 



3H- 



(10) 



The evolution equation for the peculiar gravitational po- 
tential tp is given by 



<P 



3H{1 + w)<p' - wV 2 <p = 



(11) 



whose solution in the radiation dominated era in Fourier 
space reads 



n \ 3ji(x) 



(12) 



where x = krj/y3 and ji denote spherical Bessel func- 
tions of order I. The latter expression manifests the well- 
known stagnation effect for perturbations which crossed 
the Hubble radius during the radiation era. 

We can now evaluate in Eq. Q the contribution from 
(3{rf), given by Eq. @. With the solutions Eqs. JTDJl it 
reads 



a 1 



drje ijk I 



f2d 3 <pd k tp / 



H 



7d j ipV 2 d k ip d^ipd k V 2 (p' 

12ft 2 12ft 3 
d j (p'd k \7 2 (p dV^VV 



I2ft 3 



12H 4 



(13) 



This expression can be integrated analytically exploiting 
the non-trivial result given by Kelvin's circulation theo- 
rem (see Appendix^}, which states the conservation of 
angular momentum to all orders. In fact, it can be seen 
that at second order 



1 



12ft' 



(14) 



where we simply used Eqs. (|1UI11|) . with w = 1/3, in 
Eq. (|AT0)l . assuming vanishing vorticity of the radiation 
fluid, u>i = 0, which fixes the value of our integration 
constant /3m j. One can explicitly check that the time 
derivative of Eq. (|I4|I gives the integrand in Eq. I|13|l . 
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III. POWER-SPECTRUM OF THE SEED 
MAGNETIC FIELD 

We can now compute the power-spectrum of the mag- 
netic field generated up to recombination. The correla- 
tion of the Fourier modes of the field can generally be 
written as 0] 



(^(k)^(k')) 



(27T) 



-<5 3 (k-k')[(4 



k 3 A{k)] 



kik m )S{k) 
(15) 



where k denotes the unit vector in direction of k. The 
term proportional to A(k) represents a "helical" com- 
ponent, that is a non- vanishing field component in the 
direction of the current (B • (V x B / 0)). It can be 
seen that the seed field obtained in the previous section 
has no helical component, thus A(k) = 0. The relevant 
component S(k) can be obtained from 

(B(k) ■ B*(k')> = (27r) 3 S(fc)(5 3 (k - k'). (16) 

By replacing Eq. I|14|) into Eq. (J7J), the Fourier modes of 
B can be written as 



B(M) = - 



to p (1 + z) f d 3 k 



eH 2 

2^(|k-k'|Mfc')- 

k' 2 



-k x k' 



(2tt) 3 



_/(|k-k'|M*') 



12H 



^(|k-k'|Mfc') 



(17) 



where we have neglected the last term in Eq. (J7J. 

Up to the matter-radiation equality time the gravita- 
tional potential and its derivative are given by Eq. (|12fl . 
and the correlation function of B modes can be computed 
in terms of that of <^o(k), which is a Gaussian random 
field with auto-correlation function 



(^o(k)^o(k')) = (2w) 3 P v (k)6 3 (k - k') 



(18) 



where P<p(k) is the gravitational potential power- 
spectrum, P v (k) = P 0ip k^ 3 (k/ko) ris ^ 1 and k is some 
pivot wave-number. For the normalization of P v we 
can use its relation with the power-spectrum of the co- 
moving curvature perturbation (see e.g. Ref. |l6j| ') 1Z, 
namely A^(fc) = A^(A; )(A:/A:o)™ s " 1 , where A|,(fc ) = 
(25/9)(P 0v /2tt 2 ) [13. Then, the resulting magnetic field 
is a chi-square distributed random field with power- 
spectrum 



S(k,v) 



162m 2 (l + z) 2 



e 2 H 



d 3 k' 

(27)3 *v 



P v (\k-k'\)P v (k') 



Ik x k 



4x9 



ji(xi)ji{x 2 ) 



+./j(-oi/.i.(-n 1 ( |- - X 



(19) 



o 



1e-43 



„ 1e-44 r 



% 1e-45 



1e-46 



1e-47 



0.001 




k [Mpc'] 



FIG. 1: Power-spectrum of the magnetic field at re- 
combination as a function of the comoving wavenum- 
ber. The upper and lower curves refer to the case 
with or without Silk-damping included, respectively. 



where we have defined x\ = k'fj/y/3, xi = |k — k'|f|/v3. 

The main contribution to the power-spectrum comes 
from the term in B with a time derivative of (p (second 
term in Eq. (|19l) ). As <p tends to a constant after equal- 
ity, the contribution to the generation of B during the 
period from equality to recombination is subdominant 
with respect to that up to equality time. 

We show in Figure 1 the power-spectrum of B at 
recombination, S(k), obtained from the numerical in- 
tegration of Eq. i|19[l . a s a function of the comoving 
wavenumber k. We have used the following values for the 
model parameters: A^(fco) = 2.3 x 10~ 9 , zko = 3454, 
z rec = 1088, h = 0.7 and n s = 1 [ij, EH- At small 
scales diffusion damps the fluctuations in the photon and 
baryon fluids. This effect can be taken into account by 
a factor exp(— k 2 /k%) multiplying the velocity perturba- 
tions (see e.g. [jj|)- As shown in the plots, Silk damping 
affects the produced magnetic field on scales smaller than 
about 2 comoving Mpc. 

For small k (k < 0.04 Mpc -1 ), S(k) oc k 2 , as pointed 
out in Ref. [14J, in order for the magnetic field to be 
divergenceless. For wavenumbers k > 1 Mpc -1 the mag- 
netic field spectral index approaches instead n ~ — 1, 
i.e. S{k) oc k , for scale-invariant (ns = 1) primordial 
perturbations. For a general scalar spectral index n$, 
S(k) oc k 2ns - 3 . 

The mean square value of the field on a given scale A is 
obtained by averaging over a volume of size oc A 3 0,0], 
convolving with a Gaussian window function, namely 



d 3 k 



5*(fc)exp(-A 2 fc 2 /2)- 



(20) 



For comoving scales larger than ~ 1 Mpc, up to 50 Mpc 
the rms magnetic field at recombination is well approxi- 
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mated by 



BxiVrec) « 10- 23 (A/Mpc)- 2 G. 



(21) 



After recombination the magnetic field can be considered 
to be frozen into the plasma and thus it redshifts with the 
expansion of the Universe as B(?/) = B(77 rec )(a/a rcc )~ 2 - 
Then, the average field scaled to its value today results 



B A (ry )«10- 29 (A/Mpc)^G 



(22) 



IV. CONCLUSIONS 

In this paper we have discussed a new mechanism for 
the generation of the cosmic seed magnetic field. It acts 
during the evolution of the Universe up to the epoch of 
hydrogen recombination. The underlying physical phe- 
nomenon is the vorticity of charged particles driven by 
gravitational vector modes. The latter, in turn, arise 
from the non-linear evolution of purely scalar (density) 
primordial perturbations. Hence, our mechanism is a 
generic prediction of the standard hierarchical structure 
formation scenario and does not require any ad-hoc as- 
sumption. 

Maybe the most important feature of the created mag- 
netic field discussed in this paper is that its power- 
spectrum is fully determined by the power-spectrum of 
the primordial density perturbations. As a consequence, 
a significant signal is expected over cosmological scales 
much greater than those encountered in other mecha- 
nisms acting during the early evolution of the Universe. 



Acknowledgments 

We would like to thank Klaus Dolag, Dario Grasso, 
Diego Harari and Igor Tkatchev for discussions. 



APPENDIX A: VORTICITY CONSERVATION 

We want to illustrate here how one can generate a non- 
zero rotational component of the fluid velocity at second 
order in the (dominant) radiation fluid without violating 
the conservation of vorticity, dictated by Kelvin's circu- 
lation theorem, which is non-perturbative and applies to 
a perfect fluid with equation of state p = p{p) coupled to 
gravity (see Ref. [lflj)- 

As we can see from Eq. (JSJ, we were indeed able to gen- 
erate at second order a non-zero value for /?| 7 ^ (and also 

for C| 7) ). And then, through Eq. 0, this & ^ C- p) ^ 
acts as a source for the magnetic field. The crucial point 
is that none of these quantities (neither Q nor /%) coin- 
cides at second order with the quantity that is conserved, 
according to Kelvin's theorem. 

We recall what the theorem states [l|| . First of all one 
considers a system with a single perfect fluid (in our case 



this is the radiation fluid, in the limit in which the other 
components are negligible), with four- velocity u M . Then 
one constructs the quantity 



(Al) 



where h pv = <7 M „ + u^u^ is the projection tensor into 
the rest-frame of an observer moving with 4- velocity . 
Then we split v pv it into its symmetric and antisymmetric 
parts, namely 



J p,v w lit* r ^fiu j 



(A2) 



where the antisymmetric part ui^v is called the vorticity 
tensor. The symmetric part 9^ v can be further split into 
its trace {volume expansion) and trace- free part {shear 
tensor): 



(A3) 



Finally, one can construct the vorticity vector as |23j : 
1 



-rr ' U v UJ pa . 



(A4) 



It is possible to show now, using the field equations 
[l9j . that this quantity obeys the equation 

h\{l 2 ^) = <(ZV) + l -^ Uv u p , a , (A5) 

where the overdot stands for the convective time deriva- 
tive along u M and the scalar I is defined by 



I 1 
1 = 3' 



(A6) 



thus representing the local expansion of the fluid in its 
rest frame (which coincides with the scale-factor in an 
exactly homogeneous and isotropic FRW model) . Substi- 
tuting into Eq. I|A5(I the momentum-conservation equa- 
tion one can find (see Ref. |2(j) that the following equa- 
tion holds for oj = y/uJJJIJt 1 : 



id 



3 



(P + P) au^Mhy 
= — 



(A7) 



for a perfect fluid with equation of state p = p(p), im- 
mediately leading to the result that if u is zero at some 
given initial time, then it will be always zero. 

If we deal with w as a small perturbation, the right- 
hand side of Eq. (|A7(I can be neglected (being of higher 
order due to the presence of a^"), so, we get the usual 
angular-momentum conservation law (e.g. Ref. 



lo(p + p)a 5 = constant 



(A8) 



which means that, even if we start with a small non-zero 
lo, it will decay with time at a rate which depends on the 
equation of state p{p). 
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We can show indeed that at first order w ! is related to 
the curl of the fluid velocity 



<"(!) = - 



1 

2^2' 



ijk 



(djv k + djXk) 



(A9) 



and it is well known that this vortical component is not 
generated by standard perturbation generating processes 
- such as inflation - so it can be safely set to zero. 

At second order, though, the quantity lu 1 does not coin- 
cide with the curl of the velocity £* , nor with the quantity 
(3\ but it also contains squared first-order terms |24|. 



1 



= -^[(a/3 l ) + ^pvjdkip + v 3 v' k )} 



(A10) 



This shows that in order for oj l to vanish /3 1 has to be 
generated. In other words, the momentum-conservation 
equation, whose general form for a perfect fluid at second 
order in the Poisson gauge is 



{vi + XiY + (1 - 3w)H(vi + X i) + 



w 



1 + w 



-6M 



+di<j> -\ (3tp — S)diS — (2 — 3w)cp'vi 

1 + w 



-wvfvi + -(v 2 ) } 



1 



(V 2 



j -. ■ 2 v- - 2 
leads to a conservation equation for u>. 



0, 



uj[ + 3(1 - w)Huj 1 = , 



(All) 



(A12) 



which is of course equivalent to Eq. I|A8() above. In partic- 
ular, it can be checked that the direct solution, Eq. (|13fl . 
of Eq. ijAllfl (with w = 1/3) coincides with the ex- 
pression of (3i which is found using the fact that to 1 in 
Eq. IjAlOjl is always zero. This result is given in Eq. (|14J) . 

Finally we can understand why the magnetic field is 
generated by looking at Eq. (JSJ. In fact, from there, we 



see that the magnetic field is sensitive to the vorticity 
of the charged matter components like protons. In a 
Universe with only self-gravitating pressureless matter, 
this would be conserved, so no magnetic field would be 
created. Note also that in the case of pure matter it 
is possible to show that the conserved vorticity vector 
is aM m >* = -a 2 f3^ i /2 (see Ref. 12j). Then, let us 
consider what happens to a sub-dominant matter fluid in 
a radiation dominated Universe. We may consider first, 
as an illustration, the case of a non-interacting matter 
component. Since we are in the radiation era the matter 
components yield a negligible contribution to gravity, and 
the potentials are driven by radiation. So the quantity 
p{m)i j g n0 i on g er conserved but has a source (which 
is non-zero as long as tp' is not-zero and its gradient is 
not parallel with that of (p). In fact it has to satisfy 
(the curl of) the momentum-conservation equation for a 
pressureless fluid, which is 



2 ^W (m) ) _ Q 



(A13) 



Finally, accounting for an interacting fluid (and so «W ~ 
«W and @W ~ ^W), Eq. {Hp is no longer satisfied. 
However, the left-hand side is compensated by the rota- 
tional part of a non-zero electromagnetic field, as one can 
see in Eq. @ , that we rewrite (neglecting once again the 
diffusion term and taking ~ v^ e ') as 



(m)x 



ikd j ((p'v k 0™)) 



e (a 2 i?,)' 



in 



(A14) 
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